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Abstract 

In non-integrable Hamiltonian systems with mixed phase space and discrete symmetries, sequences 
of pitchfork bifurcations of periodic orbits pave the way from integrability to chaos. In extending 
the semiclassical trace formula for the spectral density, we develop a uniform approximation for the 
combined contribution of pitchfork bifurcation pairs. For a two-dimensional double-well potential 
^T)'- and the familiar Henon-Heiles potential, we obtain very good agreement with exact quantum- 

i mechanical calculations. We also consider the integrable limit of the scenario which corresponds 

to the bifurcation of a torus from an isolated periodic orbit. For the separable version of the 
Henon-Heiles system we give an analytical uniform trace formula, which also yields the correct 



o 



, harmonic-oscillator SU(2) limit at low energies, and obtain excellent agreement with the slightly 

Q I coarse-grained quantum- mechanical density of states. 

o 

1 Introduction 



Q 
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The goal of the research area called "quantum chaos" is to relate the quantum-mechanical and clas- 
sical properties of a classically chaotic system. For autonomous Hamiltonian systems, the eigenvalue 
^ ', spectrum is, to the leading orders in h, dominated by the periodic orbits of the classical system. 

For chaotic systems, the periodic orbits are isolated in phase space and contribute individually to 
the semiclassical spectral density E] , while in integrable systems the leading contributions come 
from families of degenerate orbits forming rational tori [3 Hj. The most general case is that of 
^ I a system which is neither integrable nor ergodic but exhibits a mixed phase space consisting of 

^ ■ regular islands separated by chaotic domains. The chaotic regions increase through the destruction 

(<~^ . of rational tori when continuous symmetries are broken, and through bifurcations of periodic orbits 

OO i when the energy or another control parameter of the system is increased. Explicit semiclassical trace 

' formulae have been given for various systems with continous symmetries El El Ej ; for symmetry 

^ . breaking through the destruction of rational tori [HI El El E] i and for isolated bifurcations 

jl3l 1141 [T3] . However more complicated bifurcation scenarios which usually occur in realistic phys- 
ical systems and, in particular, bifurcation cascades I17 | ITH l I19j still constitute one of the most 
^ ■ serious problems of the semiclassical theory. 

^ , Periodic orbits contribute to the semiclassical density of states individually only as long as 

k> I they remain isolated in phase space, i.e., as long as their actions differ by large multiples of h. 

^ ■ Near bifurcations this condition is violated and the standard remedy is to determine a collective 



contribution of all periodic orbits participating in the bifurcation. In the neighbourhood of a 
bifurcation, this was achieved in jSl E] using the theory of normal forms based on the classification 
j2Ul I21j of generic bifurcations with codimension one (i.e., bifurcations occurring when one control 
parameter is varied). In all these classes, a central orbit of period n is surrounded by m > 1 
satellite orbits of period nm. The corresponding generic bifurcations are called isochronous (m = 1), 
period-doubling (m = 2), period-tripling (m = 3), period-quadrupling (m = 4), etc. The "local" 
uniform approximations developed in |B] fail at large distances from the bifurcations where the 
orbits become isolated. In |131 1141 115j "global" uniform approximations were developed, which 
interpolate between the collective contribution of the orbit cluster near a bifurcation and the sum of 
individual contributions of the isolated orbits far from it, as correctly described by the Gutzwiller 
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trace formula jlj. These global uniform approximations can, with minimal modifications, also be 
applied to nongeneric bifurcations in systems with discrete symmetries |22| . Similar global uniform 
approximations have also been derived for nongeneric bifurcations of codimension two |231 124j . Even 
though such bifurcations occur only when two control parameters meet the bifurcation conditions 
simultaneously, they are more generally of relevance because they may appear as sequences of generic 
bifurcations when one of the two parameters is fixed and only the other is varied [SHI- In other words: 
when two generic bifurcations lie so close that the orbits do not become isolated between them and 
hence the corresponding generic codimension-one global uniform approximations cannot be used, a 
description using codimension two (or higher) becomes necessary. Such a description was first given 
in |24j for codimension two along with a classification of the possible generic bifurcation sequences 
according to catastrophe theory. 

In the present paper we study a sequence of two successive isochronous pitchfork bifurcations 
of an isolated periodic orbit. This scenario which occurs in systems with discrete symmetries is not 
included in the classification of codimension-two bifurcations |23| so that at present there exists no 
semiclassical approach for it. In fact, it may constitute the beginning of a bifurcation cascade in 
which this sequence is repeated infinitely often. Such a cascade can form a geometric progression 
reminiscent of the Feigenbaum scenario [201 (although there the bifurcations are generically period 
doubling), and the new periodic orbits born at the bifurcations may exhibit self-similarity properties 
|18| I19j . Bifurcation cascades are frequently found in physical systems with discrete symmetries and 
mixed classical dynamics jl61 117| , so that a semiclassical approach to those situations would seem 
to have been required long ago. Here we develop a uniform approximation of codimension two for 
the contribution of a pair of pitchfork bifurcations to the semiclassical density of states and test 
it numerically by comparison with exact quantum-mechanical calculations. The agreement turns 
out to be very good. The degenerate limit, in which the two pitchfork bifurcations coalesce, occurs 
generically in integrable systems: there a whole family of degenerate orbits, forming a torus, is 
born from the central orbit at the bifurcation. For this case we can give analytical expressions for 
our uniform approximation, and numerical calculations for a separable system yield an excellent 
semiclassical approximation to the exact quantum- mechanical density of states. 

Our paper is organized as follows. In section 2 we present our new uniform approximation, 
whose detailed derivation is given in appendix A. The uniform approximation for the bifurcation 
of a torus from an isolated orbit in the separable limit is discussed in section 3, with its detailed 
derivation given in Appendix B. In section 4 we apply our results to a two-dimensional double- well 
potential and to the familiar Henon-Heiles system ^|, as well as to its separable version, and 
compare them to results of exact quantum calculations. An alternative derivation of our uniform 
approximation for the separable limit from EBK quantization is given in appendix C. 

2 Uniform approximation in the non-integrable case 

The density of states of an autonomous system with Hamiltonian H is given by the trace of the 
retarded Green function G{E) 

g{E) ^Y.6{E-E^) = -^Qm TrG{E) , G{E) = ^ . (1) 

n 

As usual, we split g{E) into a smooth and an oscillating part: 

g{E) = g{E) + 5g{E) . (2) 

The smooth part g{E), which semiclassically is determined by all periodic orbits of the classical 
system with zero length |28j . may either be determined by the (extended) Thomas- Fermi (TF) 
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model [ini or, where this is not analytically possible, by a numerical Strutinsky averaging of the 
quantum spectrum j291 l3Uj . The periodic orbits of finite length make up the oscillating part 5g{E). 

The semiclassical contribution to Sg{E) of any region 17 on a Poincare surface of section (PSS) 
of the phase space is given by [HI El 
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Here q' are the final coordinates and p the initial momenta on the PSS transverse to a periodic 
orbit with period T centered in the origin. The n-th iterate of the Poincare map is given by its 
generating function S {q',p,E), and the usual canonical relations hold: 
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The periodic orbits are the solutions of 

dS_ 
dq' 



dS_ 
dp 



(4) 



(5) 



which are the stationary points of the phase in If the integrals in ^ are calculated in the 
stationary-phase approximation, one obtains the individual Gutzwiller contributions of the periodic 
orbits ^ within fi: 

6g^ {E) = ^5 {E) cos - J/.^) , (6) 



(7) 



, n 2 . 

where for a two-dimensional system the amplitudes have the form 



TrM. 



The quantities S^, T^, n^, and /i^ are the action, period, repetition number, stability matrix 
and Maslov index of the orbit ^, respectively. The stationary-phase approximation yields good 
results only if the periodic orbits are isolated in phase space. Near a bifurcation this condition is 
not fulfilled so that one has to perform the integrals in © collectively over the whole periodic orbit 
cluster involved in the bifurcation. To this purpose, one inserts the normal form of the generating 
function S {q',p,E) into Q and solves the resulting integrals exactly. 

A sequence of two period-doubling bifurcations of periodic orbits is not generic because it would 
imply a jump in the stability of the central periodic orbit On the other hand, an isochronous 
bifurcation creating a new orbit with a degeneracy factor of two is equivalent to a generic period- 
doubling bifurcation The degeneracy factor two has to originate from a two- fold discrete 
symmetry of the system. Due to the behaviour of TrM^ near the bifurcation, a generic period- 
doubling bifurcation is often called a pitchfork bifurcation. The case of interest here is a sequence 
of two such nongeneric pitchfork bifurcations which can arise successively from the same central 
periodic orbit in systems with discrete symmetries such as studied in |181 119j . For this scenario we 
propose the new normal form 



S{q',P,E) = SoiE) + q'p-^{e,p' + eW')-l{p'+q" 
= So (E) + q'p — ( ei cos^ 4> + ^2 sin^ I — o,!^ ■ 



(8) 



Here q' = v2Isin(/) and p = \/2Icos<j) define the polar coordinates (/,(/>) on the PSS, with the 
central orbit sitting at / = 0. The parameters ei and €2 measure the distance to the bifurcations 
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and become zero at the first and second bifurcation, respectively. Inserting the angles = and 
(/) = 7r/2 into ((HI), one obtains the respective generic normal forms of the period-doubling bifurcations 
[Tl| corresponding to two cusp catastrophes: 



S {q' (0, 1) ,p (0, 1),E)- So {E) - q' (0, 1) p (0, /) = -^p^ _ 



(9) 



and 



S[q 



So{E) 



4" 



,/4 



(10) 

as well. The dependence of 



The stationary points of © and are the stationary points of 
the topology of (jSJ on the parameters ej is sketched in Fig. ^ 

The period-doubling bifurcations always have a real side where the central orbit as well as its 
satellite orbits are real, and a complex side where the central orbit is real but the satellite orbits 
are complex ghost orbit. We introduce a parameter cjj which is +1 on the real side and -1 on the 
complex side of the pitchfork bifurcation % with i = 1,2. Additionally, the sign of the difference 
between the actions Si of the new satellite orbits and the action 5*0 of the central orbit is indicated 
by Oi = Sign (AS'j) and A5j = Si — S^ with i = 1, 2. 






P p P 

Figure 1: Contour plots of the normal form (jSJ in dependence of the parameters for the case 
a = —1. From left to right: 62 < ei < 0, e2 < < ei and < £2 < ei. 

The uniform approximation describing the contribution of the orbit cluster involved in the 
bifurcation sequence is derived in detail in appendix A. It reads 



5giE) 
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2tt 



iir'^h- I JO 
where the functions Fi {(p) with i = 0, 1, 2 are given by 



aoFo + aiFi {4>) + a2i^2 (0)] 




C 



(11) 

(12) 
(13) 

a,Sign(e(0)). (14) 

For the evaluation of (|T2|) - dJ, any of the di can be used due to their identical values as described 
in appendix A. The coefficents oq, ai and 02 are given as solutions of the linear system of equations 
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(15) 



-2eie2 + 2e-|^ I /i /ty |— icic2 -r ^C2| 

where the amplitudes Ai with i = 0, 1,2 are given in ((Jj). C {%) and S (x) are the standard Fresnel 
functions . The index ly is related to the Maslov index /io of the central orbit by 

Z. = /X0-(cTl+CJ2)/2. (16) 

All coefficients in pi|l are expressed by the quantities which appear also in the Gutzwiller contribu- 
tions ©, which means that the uniform approximation is invariant under canonical transformations. 
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3 Uniform approximation for the separable limit 



In the degenerate case e = ei = 62 the normal form (jHJ becomes 



S{q\p,E) = So{E)+q'p-e 




2 

= 5o (^) + g'p - e/ - al2 . (17) 

Here S{q',p^E) — q'p is independent of the angle (p, which means that it refers to an integrable 
system in which the Hamiltonian depends only on the action variables but not on the angles. The 
normal form ()17|) and the corresponding bifurcation scenario has been studied in earlier works 
[HI El II '^1 132j . What we intend here is to solve the necessary integrals analytically and express all 
the coefficients by the actions and the Gutzwiller or Berry- Tabor amplitudes of the periodic orbits, 
in order to give a final formula which is easy for implementation in actual examples. 

The stationary point of the function S — q'p now corresponds to a family of periodic orbits, i.e., 
a rational torus which is created from the central orbit at the bifurcation [S]. It consists of real 
periodic orbits on one side, whereas its periodic orbits have complex coordinates on the opposite 
side of the bifurcation. To distinguish between the two sides we introduce a parameter a which 
takes the value +1 on the side where the torus is real and —1 on the side where it is complex. 

The uniform approximation for this degenerate limit ei = £2 of is derived in appendix B 
and can be given, to the leading orders h, in analytical form as 

69 (E) = ^5?e<|e^(^^--H(4=e-^^f+a 





where So and St are the actions of the central orbit and the torus, respectively, and their difference 
is denoted by 

AS = St- So (19) 

with a = Sign (AS*) = Sign(a). The amplitude ^0 corresponds to the Gutzwiller amplitude Q 
of the central periodic orbit, whereas for the torus one has to use the Berry- Tabor amplitude At 
[11IS2]- The Morse index i' appearing in p8j) is related to the Maslov index /^q of the central periodic 
orbit by 

u = fio + <7<5'- (20) 



4 Numerical results 

In order to test the above uniform approximations we apply them to two model systems: (i) a 
double-well potential which classically possesses two bifurcation cascades, one for approaching the 
saddle from below and one for approaching it from above, (ii) the familiar Henon-Heiles system as 
well as its separable version. We compare with results from exact quantum calculations and discuss 
the validity of the uniform approximations. 

4.1 A two-dimensional double- well potential 

We study the following Hamiltonian with a double-well potential 

«4(pi+p;)+i(^^-.V<'''-5-v)+^. (21) 
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Figure 2: Scaled double- well potential. Left: contour plot with the four shortest periodic orbits A 
and B evaluated at e = 0.96. Right: cut of the potential along u = 0. 



The potential in (|21() has two minima at x = and y = ibl/2\/A with energy E = 0, separated 
by a saddle at a; = y = with energy E* = 1/16A. Using scaled variables u = \/Ax and v = ^/Xy, 
the classical dynamics of the system only depends on one scaled energy variable e = E/E* = 16\E, 
with the central saddle at the height e = 1 (see figure EJ. At a scaled energy e = 9, the system 
possesses four other saddles at u = ±1 and u = ±-v/3, over which a particle can escape. At all 
energies e > 0, there exist orbits A and B that librate along and across the v axis, respectively. 
Orbit B is stable up to e = 4.778 and orbit A undergoes two bifurcation cascades, one approaching 
the saddle at e = 1 from below, and one approaching e = 1 from above. We consider here only 
energies e < 1, for which all periodic orbits appear twice corresponding to the two potential wells. 
In this region, the influence of the continuum above e = 9 can be safely neglected and the quantum 
spectrum is real and discrete to a very good approximation. We have obtained it numerically by 
diagonalisation of H21() in a finite harmonic-oscillator basis. 

The period and action of the A orbit are given analytically in terms of its two turning points, 
for V > given by 



1 r. ^ 1 



The period becomes 



and the action is 



vi = -^l-V^, V2 = -^l + V^. (e<l) (22) 
TME) = —Kiq), (23) 

V2 



Q tjr^ 2^2 
Sa{E) = ^V2 



(24) 

where E and K are the complete elliptic integrals (31j with modulus q: 



q = ^^vj-vl (25) 

The average (TF) level density of this system (including a factor 2 which accounts for the two wells) 
is given by the integral 

9a/2 \ (fn — — f?) 

grHE) = ^^ Vi^-=i^^d. (26) 

which we could not express in a simple closed form and therefore integrated numerically. 

At the energy e = 0.91232 orbit A becomes unstable, creating a stable rotational orbit R with 
Maslov index 5. At e = 0.94272 orbit A becomes stable again, creating an unstable librational orbit 
L with Maslov index 6. In figure|31the periodic orbits R and L are shown together with their complex 
"ghost" predecessors which correspond to librations in the real and imaginary parts, respectively. 
The bifurcation scenario is seen in the upper left panel of figure |1] in terms of the stability traces. 
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Figure 3: Orbits participating in a pitchfork bifurcation sequence in the double- well potential H21() . 
Upper row: Real part (left) and imaginary part (middle) of ghost orbit R at e = 0.90864, and real 
orbit R at e = 0.95 (right). Lower row: Real part (left) and imaginary part (middle) of ghost orbit 
L at e = 0.94, and real orbit L at e = 0.95 (right). 
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Figure 4: Pr operties of the periodic orbits A, R and L near their bifurcations in the double-well 
potential (PT|) . plotted versus the scaled energy e. Top left: stability traces; middle left: action 
differences; bottom left: periods; and right: Gutzwiller amplitudes (cf. text). The dashed portions 
of all curves correspond to the complex pre-bifurcation ghost orbits. 

In figure |3 we also show the action differences and periods of the three orbits, as well as their 
Gutzwiller amplitudes, plotted versus the scaled energy e. As shown analytically in [Til, the asymp- 
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Figure 5: Oscillating part of density of states in the double- well potential (PT|) . Solid line: exact 
quantum result obtained with A = 0.0008. Dashed line: uniform approximation including isolated 
contribution of orbit B. Dotted line: sum of Gutzwiller contributions of isolated orbits, diverging 
at the two lowest bifurcations of the A orbit. (The other bifurcations, lying at e > 0.9998, cannot 
be seen at this resolution.) Coarse-graining by Gaussian convolution with energy width 7 = 0.5. 



totic divergences of the amplitudes of the central orbit (here A) and the satellite orbits (here R and 
L) must differ by a factor \/2 ; this factor has been included in the right panel of figure in order 
to confirm this fact numerically. 

Using these numerical results we now evaluate the uniform approximation (jllj) for the joint 
contribution of the orbits A, R and L. The B orbits are included in the standard Gutzwiller approx- 
imation, since they stay isolated at all energies and do not interfere with the other orbits. In figure 
inithe result is shown together with the result of an exact quantum-mechanical diagonalization done 
for A = 0.0008. One can recognize that the uniform approximation tremendously improves over 
the diverging standard Gutzwiller approximation (dotted line), leading to a perfect agreement with 
quantum mechanics up to the saddle at e = 1. Here, as well as in all following comparisons with 
quantum mechanics, we have coarse-grained the density of states by convolution with a Gaussian 
over an energy interval 7. In the semiclassical trace formulae this leads i29j to the inclusion of an 
exponential factor exp{ — (7Tg/2?i)'^} in the Gutzwiller amplitude of each periodic orbit ^, where 
is its period, in regions far enough from the bifurcations for the orbits to be isolated. Note that 
in the regions between the two bifurcations, the Gutzwiller approximation is not valid, so that our 
codimension-two uniform approximation is indispensible. 

4.2 The Henon-Heiles system 

The system of Henon and Heiles is given by the Hamiltonian Wl\ 



H = \ (vl + P^) + ^ (f + y')+x Uv - ^ y^) . (27) 



When the scaled variables u = Xx and v = Xy are introduced, the scaled total energy in units of 
the saddle-point energy E* = 1/6A^ becomes 



e = E/E* = 6 



- (li^ + + y (u, v) = 3 (u^ +3 (n^ -Fu^) -^6^2 - 2t>^ (28) 



In the left part of figure El we show the equipotential lines of the potential part of (|28() in the 
(n,u)-plane together with the three shortest periodic orbits A, B and C, evaluated at the scaled 
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energy e = 1. Along the tree mirror axes (dashed hnes) the potential is a cubic parabola as shown 
along n = in the right part of figure IHl For an arbitrary energy e < 1, the turning points of the A 
orbit are determined as the solutions of the equation 

e = 3-^2 + 2v^ (29) 

and given by 

t;^ = l/2-cos(7r/3-</>/3) , W2 = 1/2 - cos (vr/S + 0/3) , = 1/2 - cos ((/>/3) , (30) 

with 

cos ((/>) = 1 - 2e . (e < 1) (31) 

As was shown by Henon and Heiles, the classical dynamics is quasi-regular up to energies of about 
e = 2/3 and then becomes increasingly chaotic |22j. The v motion of the A orbit with the scaled 
energy e„ is given by jTHj 

v^A^\t) =vi + {V2 - vi) sn2 {s, q) (32) 

in terms of the Jacobi elliptic function 31 sn(s,g) which depends on the argument s and the 
modulus g, given by 



s = t\l (w3 — f i)/6 and Q = \ — (33) 

The turning points Vi have to be evaluated according to the equations (|5n|l with e = Cy. 




Figure 6: The Henon-Heiles potential. Left: Equipotential contour lines in scaled energy 
units e in the plane of scaled variables u, v. The dashed lines are the symmetry axes. The 
three shortest periodic orbits A, B, and C (evaluated at the energy e = 1) are shown by 
the heavy solid lines. Right: Cut of the scaled potential along u = 0. 

The periodic orbits of the system have been investigated and classified by Churchill et al. j33j 
as well as Davies et al. |31]. Up to energies of e w 0.97 there exist only three types of periodic 
orbits with periods T of the order of 27r: the librations A and B, and the rotation C. Due to the 

symmetry of the potential, orbits A and B occur in three orientations connected by rotations 
about 27r/3 and 47r/3 in the {u,v) plane. Orbit C has a degeneracy of 2 because of the time reversal 
symmetry which corresponds to two different orbits with opposite senses of rotation. The orbit B is 
unstable for all energies and the orbit C stays stable for energies below e = 0.8922 where it becomes 
unstable due to a generic period-doubling bifurcation. The bifurcation cascades of the A orbit and 
the orbits generated by them have been studied in detail in jl81ll9j : we adapt the names of the orbits 
given in these references, whereby the subscripts of the orbit names denote their Maslov indices. 
The A orbit undergoes its first isochronous pitchfork bifurcation at an energy ei = 0.969309 and the 
second one at 62 = 0.986709. At the first bifurcation it creates a stable rotational orbit R5 which is 
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doubly degenerate due to its two possible senses of rotation. At the second bifurcation, it creates 
an unstable librational orbit Lg which is doubly degenerate due to the reflection symmetry of the 
potential at the v axis. This scenario repeats itself at higher energies, whereby the pairs of orbits R7 
and Lg, Rg and Lio, etc, are born. The rotational or librational character of these orbits is indicated 
by the letters R and L, respectively. In figureI7|the orbits R5 and Lg are plotted together with their 
pre-bifurcation complex ghost orbits. All orbits, including A, gain one more degeneracy factor 3 
due to the three-fold discrete rotational symmetry of the potential, so that the overall degeneracy 
factors of A, R and L orbits are 3, 6 and 6 respectively. 




Figure 7: Orbits born in the first pitchfork bifurcation sequence in the Henon-Heiles potential. 
Upper row: Real part (left) and imaginary part (middle) of ghost orbit R5 at e = 0.9690, and real 
orbit R5 at e = 0.9798 (right). Lower row: Real part (left) and imaginary part (middle) of ghost 
orbit Le at e = 0.9864, and real orbit Lg at e = 0.9870 (right). 




u u u 

Figure 8: Poincare surfaces of section (PSS) of the scaled Henon-Heiles Hamiltonian (|28j) . taken 
for u = 0. Left: e = 0.969; middle: e = 0.982; right: e = 0.989. 

In figurelBl a part of the PSS for f = is plotted for energies before the first pitchfork bifurcation 
(left), between the first and second bifurcation (middle), as well as after the second bifurcation 
(right). The topology in the vicinity of the bifurcation sequence is correctly described by the 
normal form Q, as can be seen by a comparison with figure ^ 

In evaluating the uniform approximation, one can exploit the fact that the actions and the 

periods of the orbit A can be calculated analytically. The action is given by 

/""2 2 I 

Sa{E) = 2J Ve - 3v^ + 2v^ dv = -^^6 {v^ - vi) [E (q) + cK (q)] , (34) 
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where the modulus q of the complete elliptic integrals is given in (|33|) . The constant c is given by 



2 

c= --{vz-V2){2vi-V2-vi) (35) 
y 

in terms of the turning points Vi {i = 1,2,3) given in (|3fl|) . The period is obtained as 

T. iE) = ^ = 2^3 r = K (,) . (36) 

In figure ini the quantities needed to evaluate the uniform approximation (|llj) of the density of 
states are shown as a function of the scaled energy e. One can see that the stability trace TrM4 of 
the A orbit takes on the values +2 at the bifurcation energies. The stability traces of the orbits 
and Le are also plotted; they stay real even for energies e < ei and e < 62, respectively, where the 
two satellites are complex ghost orbits (which their properties shown by dashed lines in figure 

In jnS] , it was shown that the coarse-grained quantum- mechanical density of states of the Henon- 
Heiles potential (obtained with a Gaussian smoothing width 7 = 0.25) can be rather accurately 
approximated semiclassically, using just the isolated orbits A, B and C and their second repetitions, 
for energies far enough from the harmonic-oscillator limit e = 0. In 12 ,, a uniform approximation for 
the symmetry breaking at e = was developed which continuously interpolates from the harmonic- 
oscillator limit, given in H50() below, to the region where the Gutzwiller trace formula for the isolated 
orbits is valid. However, the bifurcations of the A orbit have not been treated uniformly in the 
references so that the accuracy of the results decreased near the saddle at e = 1. In ^S] 

the classical bifurcation cascade in the Henon-Heiles potential was discussed, in which the sequence 
of two successive pitchfork bifurcations repeats itself infinitely often. 
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Presently we test our uniform approximation to the density of states against the quantum- 
mechanical result obtained for A = 0.03. The quantum spectrum was, as in obtained by 
diagonalisation of ()27|) in a finite harmonic-oscillator basis - thus neglecting the effects of quantum 
tunnelling through the barrier. Both quantum and semiclassical results were coarse-grained with a 
Gaussian width of 7 = 0.4; for this resolution the inclusion of the second repetitions of all periodic 
orbits in (|1H) was necessary (cf. ^21)- For the pitchfork bifurcation of the second repetition of 
the orbit C at e = 0.892, where a double-loop orbit D is created jSS], we used the codimension-one 
uniform approximation of 15. The upper part of figure fTOl shows the entire energy region < e < 1, 
whereas the lower part shows the zoomed region 0.88 < e < 1. The solid lines give the quantum- 
mechanical result, and the dashed lines the results obtained with our uniform approximation Ulljl 
for the first two pitchfork bifurcations of the A orbit. In the region e < 0.5, we have included the 
uniform approximation for the symmetry breaking, developed in , in order to obtain the correct 
harmonic oscillator limit for e — > 0. The dotted line in the lower part of the figure corresponds 
to the sum of the isolated periodic orbits according to the standard Gutzwiller trace formula J]. 
Here the divergences due to the lowest bifurcations of the A and C orbits are clearly visible. The 
uniform result Ullj) . however, exhibits no divergences and its agreement with the quantum result 
is very satisfactory. The discrepancy arising at e ^ 0.992 can be attributed to the influence of the 
continuum that starts at e = 1 which was not taken properly into account in our quantum result. In 
fact, the rather excessive maximum appearing in the latter around e ^ 0.994 makes us believe that 
the latter is erroneous, rather than our semiclassical result. Note that the uniform approximation 
properly yields the asymptotic Gutzwiller result on either side of the double-pitchfork bifurcation. 
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Figure 10: Oscillating part of density of states in the Henon-Heiles potential. Solid lines: quantum- 
mechanical results obtained for A = 0.03. Dotted lines: sum of Gutzwiller contributions (jH)) of all 
isolated orbits. Dashed lines: codimension-two uniform approximation (|11|) for the orbits A, R5 and 
Lg, including orbits C and D in the codimension-one uniform approximation of 14 and the isolated 
B orbit. Coarse-graining with Gaussian width 7 = 0.4. 

In the energy region e > 1 above the barrier, where the spectrum of the Henon-Heiles Hamilto- 
nian (|27|1 is continuous, the oscillating part of the density of states is determined by the resonances 
in the continuum. In order to test the semiclassical periodic orbit theory in this domain, it becomes 
necessary to calculate both the positions and widths of the resonances. It will then be an interesting 
question to study which periodic orbits are important in the continuum region. Work along these 
lines is in progress jHEl- Although the continuum region is also classically unbounded, all the R 
and L orbits bifurcating from the A orbit (which itself ceases to exist above e = 1), as well as the 
D orbit bifurcating from C, continue to exist and are bounded at all energies e > 1 |181 134j . In 
addition, three new orbits librating across the saddles exist in this region |18 | I33 j : since they have 
the shortest periods they are expected to play a leading role in the coarse-grained density of states. 
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4.3 The separable Henon-Heiles system 

The Henon-Heiles system permits chaotic motion because of the nonseparable term x^y in (j^Tj). 



Omitting this term one obtains a system which is separable in x and y and hence integrable: 

A , 



■y 



(37) 



Again using scaled variables u 
energy E* reads 



Ax and v = Xy the scaled energy e in units of the saddle point 



e = E/E* = 6 



i (^t^ + i;^) +V{u,v) 



3 (u"^ + v"^ 



2 I 2 
U + V 



2v' 



(38) 



Figure ITTI shows a contour plot of the potential part of ()38() in the {u,v) plane together with the 
two shortest periodic orbits A and B calculated at an energy e = 1. The two orbits are librations 
along the u and v axes. The potential along the v axis is the same as that in the right part of figure 
El while the potential along the u axis is harmonic. 



> 
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Figure 11: Equipotential lines in the {u,v) plane for the separable 
version of the Henon-Heiles potential. The heavy solid lines show the 
two shortest periodic orbits A and B evaluated at e = 1. 



Th actions and periods of the A orbit are given by (|34)1 and (|36j) . respectively. The trace of its 
stability matrix is given analytically by |18j 

TvMa {E) =2 cos {Ta{E)). (39) 

The u motion of the B orbit is harmonic 

UB{t) = ^ sin {t + ^) , (40) 

where Cu is the conserved scaled energy in the u direction and the phase (p is arbitrary. The action 
and period of the primitive B orbit are 

Sb (E) = 2ttE , Tb (E) = 2tt . (41) 

The trace of the stability matrix of the B orbit has the constant value TrM^ = +2, which is 
consistent with its appearing as a torus in the asymptotic analysis given in appendix C. 
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The ky-th repetition of orbit A bifurcates whenever the condition 

TrM^" = 2 cos (k^TA) = +2 (42) 

is obeyed, which is equivalent to the resonance condition (cf. appendix C.2) at the bifurcation 
energies -E^/ 

kvTA{Ehif) = 2iTku = kuTs ■ (43) 

Thus, the bifurcations of the A orbit create the rational tori corresponding to the ky : k^ resonances. 
The new tori form families of degenerate periodic orbits that are related by the U{\) symmetry due 
to the freedom in choosing the phase (/> S [0, 2tt) in their u motion 

UT{E) = ,J^^^s\n{t + (44) 

where e^if are the scaled bifurcation energies, while their v motion is "frozen" and identical to that 
of the A orbit given in (|,'-{2|) at the corresponding bifurcation energy: 

VT{t) = v^T''\t) . (45) 

The actions of the tori become 

St{E) = k, Sa (Ebif) + k^27T {E - Euf) , (46) 
so that their periods stay constant at 

Tt = ku2^ = kuTB . (47) 
Like for all degenerate orbit families, their stability trace is constant: 

TrMr = +2 . (48) 




Figure 12: Same as figure |1] for the separable Henon-Heiles system H37() for the bifurcation of the 
fc„ : A;^ = 5 : 3 resonance at energy e = 0.987655. The central A orbit is labeled by "0", the 
bifurcated 5:3 torus by "1". 
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We first apply our uniform approximation to tlie single isolated bifurcation with ku : = 5 : 3 
which happens at e = 0.987655. In figurelT^we show the action difference Si — Sq = SxiE) — SAiE), 
the periods Tq = 3Ta{E) and Ti = Tt = IOtt, the traces of the stability matrix, as well as the 
Gutzwiller and Berry- Tabor amplitudes of the isolated A orbit and the 5:3 torus, respectively. This 
figure should be compared with figure |H1 in which the corresponding quantities are shown for the 
non-integrable Henon-Heiles potential. Here the two bifurcations coincide, and instead of the two 
isolated orbits R5 and Lg created at the two bifurcations there, we have here only one torus whose 
stability trace has the constant value +2. 

These quantities are now used to evaluate the uniform approximation for the integrable case, 
given in (|18j) . The result is shown in figure El by the dashed line. It is compared to the exact 
quantum-mechanical curve (solid line) obtained for A = 0.04, as well as to the result of including 
independently the Berry- Tabor contribution of the torus and the Gutzwiller contribution of the 
isolated A orbit which diverges at the bifurcation (dotted line). All results have been coarse- 
grained by convolution with a Gaussian with smoothing parameter 7 = 0.1. We see that the 
uniform approximation reproduces the quantum result very accurately. 
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Figure 13: Oscillating part of level density for the separable Henon-Heiles system ()37(1 near the 5:3 
resonance, coarse-grained with a Gaussian width 7 = 0.1. Solid line: quantum-mechanical result 
obtained with A = 0.04; dotted line: sum of Berry- Tabor contributio n of 5:3 torus and Gutzwiller 
contribution of isolated A orbit; dashed line: uniform approximation (|18)1 . 

So far, we have discussed and tested our uniform approximations for a double-pitchfork sequence, 
based on the normal form Q, and its separable limit. In appendix C, we give an alternative 
derivation of the uniform approximation for the separable limit, starting from the EBK quantization 
and exploiting the convolution property of the density of states for separable systems. There we 
do not require any normal form, but we start from a one-dimensional integral (jlOHj) for the density 
of states which by construction is uniform in the sense that it does not diverge at any energy. 
By expanding the amplitude and phase functions of the integrand around the bifurcation energies 
Ebif up to first and second order, respectively, we arrive at approximate integrals which precisely 
correspond to those obtained from the normal form Q, and which can be reexpressed in terms of 
the Gutzwiller amplitude of the isolated A orbit and the Berry- Tabor amplitudes of the rational 
tori. Furthermore, the starting point (|1()8)) allows us also to include the limit e — > 0, in which the 
amplitude of the isolated A orbit also diverges, in a uniform way. 

Since all amplitudes, actions and periods of the isolated A orbit and the tori bifurcating from 
it can be given analytically for the IHH potential, it poses no problem to sum over the repetitions 
of the A orbit and all the tori bifurcating from them. As shown in detail in the appendix C, 
this leads to the following "grand" uniform approximation which is valid and finite also in the 
harmonic-oscillator limit e — > 0: 
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CO oo I / \ I hj \ 

SgumiE) = E E (-l)'"+'^<| ( (E) - - ^ M^ujj cos 



+ 5g^^'\E) + 6£'\E). (49) 

All quantities appearing above are given analytically in equations l|112|) . pil-^j) . (|12()|) and (|12H) of 
the appendix C. 

The first term in ()49() yields, upon summation over all ku and and adding the term 6ga^^\E) 
in the last line, precisely the Gutzwiller trace formula (|116j) of the isolated A orbit which diverges at 
the bifurcations and at -E = 0. The second term in the first line is a counter term from the tori that 
cancels all divergences of the Gutzwiller amplitudes. The second line of (|49|) yields the Berry- Tabor 
trace formula (jlllj) far away from the bifurcations; near the bifurcations it contains the Stokes 
factor that interpolates between the Berry- Tabor amplitudes above and zero below the bifurcations, 
yielding exactly half the Berry- Tabor amplitudes at the bifurcations. The two contributions in the 
last line of (|49|) are small boundary terms, given in H114|) and H117|) of the appendix C, which are 
numerically insignificant but have been included in order to be consistent up to order in the 
amplitudes. 

In the limit e — > 0, where we can neglect all bifurcations, only the diagonal terms with ky^ — 
ky = k contribute. The trace formula ()49() then leads uniformly to the correct SU(2) harmonic 
oscillator limit whose trace formula is given in equation (|in6|) of the appendix C: 

SgumiE) Sgt%E) = ^J2^cosi^-27rEj (for e ^ 0) . (50) 

(The same limit was obtained in a uniform approximation for the full non-integrable Henon-Heiles 
potential in ^5] neglecting, however, the bifurcations.) 

In the figures im and n^l we compare the results obtained from the grand uniform approximation 
(|49|) with those of quantum- mechanical calculations for system (|38() with A = 0.04 (with saddle 
energy E* = 104.666 corresponding to e = 1), both coarse-grained by a Gaussian convolution with 
an energy range 7 = 0.1, including repetition numbers up to jfcul, < 8 into the semiclassical trace 
formula (|49|) . Figure El shows the lowest energy range which exhibits for e ^ 0.1 the harmonic- 
oscillator limit ()50() where the amplitude of ^^'(e) is linear in e. 




Figure 14: Oscillating part of level density of the separable Henon-Heiles system (|37|) . coarse-grained 
with 7 = 0.1. Solid lines: quantum-mechanical result for A = 0.04. Dashed lines: semiclassical 
results with ku,ky < 8. 
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Figure 15: Oscillating part of level density of the separable Henon-Heiles system H37() . coarse- 
grained with 7 = 0.1. Solid lines: quantum-mechanical result. Dashed lines: semiclassical results 
with ku,k^ < 8. Top: Berry- Tabor result for the tori. Centre: Sum of Berry- Tabor res ult for the 
tori plus Gutzwiller result for the isolated A orbit. Bottom: uniform approximation (|49|). 

In the top panel of figure ITKl we compare the quantum result to the standard Berry- Tabor trace 
formula, given in Hlll() of appendix C, which takes into account only the tori with semiclassical 
amplitudes proportional to h~^^'^. In the center panel, we have added to them the A orbit contri- 
bution described by the Gutzwiller trace formula, given in (|115() of appendix C, with amplitudes 
proportional to . The latter is seen to diverge at all bifurcations corresponding to resonances 
with ku : > 5 : 4:. Between the bifurcations, the result is clearly improved by adding the A orbit 
contribution and comes very close to the quantum result. In the bottom panel, finally, we show the 
grand uniform approximation (|49j) which reproduces the quantum result perfectly throughout the 
whole energy region. The bifurcation corresponding to the resonances with ku '■ k^ = 2 : 1 happens 
at the scaled energy e = 0.998491; all bifurcations with ku : k^ > 2 : 1 happen thus in the top 
0.15 percent of the energy scale very near the barrier. In this region, the bifurcations are lying so 
densely that their independent summation in (|49j) is strictly not justified. However, at the present 
resolution of the spectral density this does not appear to affect our numerical result. On the other 
hand, the good agreement which we find in figure El at all lower energies demonstrates that our 
grand uniform approximation H49|) successfully sums all partial bifurcation cascades of the A orbit 
limited by the repetition numbers 2 < ky,ku < 8. 

We should stress that, like for the non-integrable Henon-Heiles potential, the quantum spectrum 
was obtained here by diagonalisation in a finite harmonic-oscillator basis. The persistence of our 
good agreement up to e ~ 1 therefore suggests that the barrier tunnelling effects are negligible - at 
least within the resolution given here by the coarse-graining width 7 = 0.1. 
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5 Summary, Conclusions and Outlook 



We have derived a co dimension- two uniform approximation for the joint contribution of the periodic 
orbits involved in a double pitchfork bifurcation sequence by constructing a suitable normal form. 
This bifurcation scenario only occurs in systems with discrete symmetries and cannot be treated 
by the codimension-one uniform approximations developed in ^1] due to the vicinity of the two 
pitchfork bifurcations. Furthermore it does not belong to the unfoldings classified in [2^, so that a 
new approach became inevitable. We have also studied the limit where both pitchfork bifurcations 
coincide, resulting in the bifurcation of a torus from an isolated orbit such as it happens in integrable 
systems. For separable potentials, the same uniform approximation could be rederived from an EBK 
trace formula that accounts both for an isolated orbit and for the tori bifurcating from it. 

Our uniform approximation was tested numerically for two well-known systems with mixed 
classical dynamics: a double-well potential and the familiar Henon-Heiles system. In both cases the 
uniform approximation was shown to reach the asymptotic Gutzwiller approximation on either side 
of the double-pitchfork bifurcation, while yielding finite amplitudes throughout the whole energy 
region. The agreement of the semiclassical and quantum-mechanical coarse-grained level densities 
was found to be excellent. 

Our uniform approximation is only valid as long as the considered pair of pitchfork bifurcations 
is isolated from other bifurcations. In the examples studied here, this is the case for the lowest pair 
of bifurcations of the isolated A orbit. Since this orbit undergoes an infinite bifurcation cascade 
cumulating at the barrier energy e = 1, our approximation will eventually fail for higher double- 
pitchfork bifurcations; the precise energy where this happens depends on the value of the nonlinearity 
parameter A. However, for the coarse-grained shell structure obtained with a limited resolution 
(given by a sufficiently large Gaussian width 7), the higher bifurcations become less important and 
the corresponding divergences in the level density cannot be resolved. 

For the separable limit of the Henon-Heiles system we have obtained analytical expressions for 
the uniform approximation. This allowed us to sum over a large part of the bifurcation cascade 
corresponding to rational : ky tori with ku,ky < 8. The resulting grand uniform approximation 
(|49|) for the semiclassical density of states, which also correctly describes the SU(2) symmetry- 
restoring limit for e ^ 0, leads to an excellent agreement with the slightly coarse-grained quantum- 
mechanical density of states even up to the barrier energy e = 1. 

An extension of the semiclassical analysis of the density of states to the energy region above the 
barrier (e > 1), including a rigorous quantum- mechanical determination of the widths and energy 
shifts due to barrier tunnelling in the quasi-bound region of the spectrum and of the resonances in 
the continuum region, is in progress |36j . 
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Appendix A. Derivation of the uniform approximation 



In this section we describe the procedure leading to the uniform approximation following the 
ideas outlined in ^24^. The semiclassical approximation of the density of states is given by 



5g{E)'. 
with a phase function 

and an amplitude function 



5Re / dq' dp (q , p) ex.p 



-<^{q,p) -t-i^ 



^ {q',p) = S {q,p) - q'p 

1/2 



■^{q,p) 



ldS_ 

n dE 



dq'dp 



(51) 
(52) 
(53) 



The integration is done over any region Vt of the PSS and v corresponds to the Morse-Index. We 
can set n = 2 in H53() because the nongeneric pitchfork bifurcations are equivalent to generic period- 
doubling bifurcations 22 . Using the normal form Q of the generating function the phase function 
^ can be expressed in canonical polar coordinates / and (j) as 



$ {q' ((/., I) ,p{<P,I)) = So- (ei cos^ cf> + e2 sin^ a/^ 



with 

p = \p2J cos , g' = \p2J sin . 
The periodic solutions in Q correspond to the stationary points of ^ at 



yielding 



sm 



(2</>) = (ei cos^ (l) + e2 sin^ 0) + 2a/ = . 



(54) 
(55) 

(56) 
(57) 



There are four solutions of (|57j) : two with cos {2(j)) = 1, corresponding to a satellite orbit which is 
labeled by 1 in the following, and two with cos {2(j)) = — 1 corresponding to a satellite orbit labeled 
by 2. At the stationary points the values of / are 



/• = - — 

2a ^ 



(58) 



where i = 1,2. With ctj = Sign(Ii) the satellite orbit i is real if cjj = -|-1 and represents a ghost 
solution with complex coordinates q and p if cjj = —1. The phase function (|54|) evaluated at the 
stationary points H58|) corresponds to the actions of the two satellite orbits 



Si — So + . 

4a 



The periods are Tq = dSo/dE and 



To + 



ei dei 



2a dE 

The traces of the stability matrix can be calculated from 



TrM 



/ d'^s \ 

I dpdq' J 



-1 



1 (d^sV d^sd^s 

\ dpdq' J dp^ dq'^ 



(59) 
(60) 

(61) 
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evaluated at the stationary points ^lElE]- One obtains 

TrMo = 2 - eie2 , TrMi = 2 + 2eie2 - 2e? , TrMs 



2 + 2ei€2 - 2ei . 



(62) 



For all orbits the actions ()59() . periods (jbUf) and stabilities ()52() are real quantities even though the 
orbits themselves can be complex. This characteristic of period-doubling bifurcations which is due 
to a Stokes transition was already mentioned in jl4j . 

The Maslov indices /Xj of the periodic orbits are related to the Morse index u appearing in (|51|) by 



^i'i = v+-{nn- Up) , 
where and Up are the number of negative and positive eigenvalues of the matrix 



^2$ Q2q> 
9<j'2 dq'dp 



\ dq'dp 9p2 / 

evaluated at the stationary points. They follow as 

liQ = + (Sign (ei) + Sign (e2))/2, 

111 = + (Sign (e2 - ei) - Sign (ei)) /2, 

112 = + (Sign(ei - €2) - Sign(e2)) /2. 

For the amplitude function the following ansatz was found to be sufficient: 

*((/>,/) =ao + ai/ + a2^^- 
Equation (|51j) . expressed in canonical coordinates, now takes on the following form 

1 



(63) 



(64) 



(65) 
(66) 
(67) 



5g{E) 



4,2,-2 exp 





vr 


r2TT 






— l—V 


/ d4 


./ 


2 


Jo 


Jo 



X exp 



ei cos^ <A + £2 sin^ ' 



(69) 



The parameters measure the distance to the bifurcation i. They are given by the actions (|59|) of 
the new born orbits as 

a = -2aiai^\ASi\ . (70) 

where we have set 

a = ai = Sign (ASi) . (71) 

In order to achieve a uniform approximation one evaluates (|69|) in stationary-phase approximation 
at the stationary points which yields 



55^"^) {E) 



X cos 



1 (ao + aili + a2li 



h 



ei cos^ (f)i + 62 sin^ > 



Ji - ah; 



(72) 



One can now determine the coefficients ao, ai and a2 by identifying the Gutzwiller Amplitudes A: 
with 

1 (ao + o^ili + a2/j^) 



Ai 



where i = 1,2,3. Defining 



e (0) = ei cos^ </> + £2 sin^ 1 



(73) 



(74) 
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the integrals with respect to / in (|69j) can be calculated analytically using 

Jo \ de{<t))J Jo 

They yield, for n = 0, 1 and 2: 



(75) 




TTh I 1 



c 



2Trh 



idjS 



2nn 



F,ic^) = -^[zh + ei<l>)Fom , F2{<P) = -^ 



1 



(76) 



(77) 



with a = — CTjSign (e (0)). The remaining (p integral over the interval [0, 2it] can easily be calculated 
numerically. 

Appendix B. Derivation of the uniform approximation 
for the separable limit 



In the case ei = €2 ^ e the phase function simplifies to 

<f{q' {(t>,I),p{(f>,I))=So-eI-al' 



(78) 



and becomes independent of (p, corresponding to an integrable system. The stationary point of ^ 
corresponds to a torus with the radial coordinate 



It 



2a 



With the definition 



a = Sign (It) , 

the torus is real if o" = +1, while it is imaginary if cr = — 1 which can be understood from 



p2 + q'^ 



The action of the torus becomes 



and for the period one obtains 



St = Sq + 
Tt = To + 



4o' 
e de 



2a dE 



(79) 
(80) 

(81) 

(82) 
(83) 



Using (|61|) one finds that TrM = +2 which is characteristic of an orbit family. The amplitude 
function can be derived from H68|l . resulting in 



* (/) = (ao + m , 



with 



(84) 
(85) 
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using ao, ai and from (|68j) . This can be seen by the fohowing integration by parts 

~ — n dl Iexp( = — n dl I exp ( (-e - 2al + 2al) 
2a Jo \n J 2a Jo \n J 



dl /^exp(-$) +0{h), 



1 r°° ( i 

— / dl lexp 
2a Jo \ri 



2al 



(86) 



where the integral that was neglected in the last step is of relative order fi. Thus we obtain, to 
leading order in Ti, the result (|84j) . 

The functions in H74|) no longer depend on (j) so that the integration over (p can be performed 
giving a factor of 2tt. The remaining expression for 5g {E) then has the form 



5g{E) 



5Re exp 



I „ .vr 

— Do — l — V 

h 2 



/ (ao + exp --[el + ai 
Jo V Ti^ 



(87) 



Exactly the same formula can be found in ^3] in relation with a special case of a generic period- 
quadrupling bifurcation. Using the integrals (|76|) . i[77jl. one arrives at 



1 /3 /5o 
— — cos — 
vrn 2a V 



, A 1 



X 5Re < e 



c 




(3e\ 
2a j 



iaS 



2ttTi \a\ 



It remains to express all the parameters by the quantities that enter into the asymptotic contribu- 
tions of the torus and the central orbit. A stationary phase approximation of (|88|) would deliver the 
contribution of the stationary point corresponding to the torus only. In order to obtain the contri- 
bution of the central periodic orbit at I = one has to include also the end-point corrections to the 
stationary phase approximation (cf. appendix C.2). This amounts to an asymptotic expansion of 
the Fresnel functions for large arguments a; S> 1 (cf. [HJ)- Keeping their two leading terms 



leads to 



5g{E) 



+ 



C{x) 
Six) 



I ao f So 
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^ + sin (7rxV2) , 
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- — cos ynx /2 



(89) 
(90) 



vr 



u + Sign (e)] 
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a V 2a 



vr / a 
- + - 
2 V 2 



(91) 



Two asymptotic contributions can be recognized. One is of the order h~^, corresponding to the 
central periodic orbit, and one is of the order h^'^^'^ which is the torus contribution. Asymptotically 
one obtains a torus contribution only on the real side of the bifurcation a = +1 whereas on the 
complex side a = —1 the torus contribution asymptotically vanishes even though the torus amplitude 
itself must not necessarily go to zero. The fact that it still gives no contribution asymptotically is 
due to a Stokes transition of the torus. 

Expressing now ()88|) with the Gutzwiller amplitude Aa of the central orbit 



Ao = Aa 



1 ap 

nh lei 



(92) 
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and the Berry- Tabor amplitude of the torus 



^T = ^,^(ao-^), (93) 



Trn^/^y \a\ V ' 2a, 

and setting AS = St — Sq as weh as = fiQ — Sign (e) finahy yields the uniform approximation p8() . 

Appendix C. Alternative derivation of the uniform ap- 
proximation for the separable limit from EBK quantiza- 
tion 

In this appendix we give an alternative derivation of the uniform approximation for the separable 
limit, starting from EBK (or WKB) quantization in one dimension and using the fact that the two- 
dimensional density of states can be obtained by a convolution of the two one-dimensional densities 
of state. We present the general formulae, which to our knowledge have not been given before in 
the literature, in the first subsection and derive from it the known exact trace formula for harmonic 
oscillators. In the second subsection we specialize to the integrable Henon-Heiles system and present 
a uniform trace formula which sums over all bifurcations and leads to the correct harmonic-oscillator 
limit for e — > 0. 

1. Semiclassical trace formula for separable Hamiltonians 

For a separable Hamiltonian in two dimensions (u, v) 

H = Hu{u,pu)+Hy{v,py) = Eu + E^ (94) 
the Schrodinger equation separates 

^nm{u,v) = 4)n{u)^m{v), Enm = + £m (95) 

and the exact quantum density of states can be written as a convolution integral over the two level 
densities of the one-dimensional systems: 

g{E) = r gu{E - E') g,{E') dE' , (96) 

JO 

where 

g^{E) = Y,5{E - en) , g,{E) = Y,6{E - e^) , (97) 

n m 

and we have assumed em,en > 0. We now use EBK quantization (i = u, v; Uu = n, = m): 

Hi = H{Ii), I^ = l-Si = ^j>pidqi = h{ni + 1/2). (98) 

and Poisson summation (cf. [22] ) to obtain the following semiclassical trace formula for each of 
the one-dimensional level densities 



^'^B)=^§ T. (-1)'— 



(99) 
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which is identical to the Gutzwiller trace formula |1| for a one-dimensional system and yields the 
corresponding EBK (WKB) spectrum. Using H96j) . we thus get the two-dimensional trace formula 
for the separable Hamiltonian H94() 



giE) 



{2Trh) 



1 °o nE 

J2 ^.ifu+k.J T^{E - E')%{E') 



-j^T} — OO 



X COS 



COS 



Su{E-E') 



S^E') 



dE' . 



(100) 



If the convolution integral is done exactly, this trace formula yields the spectrum Enm in (US)) in 
the EBK approximation. The contribution from ku = = yields the average Thomas-Fermi 
(TF) level density which becomes a simple convolution integral over the primitive periods of the 
two one-dimensional motions: 



g{E) = grriE) = T^{E - E') %{E') dE' , 



(101) 



The semiclassical trace formula (jlOOf) . which contains the smooth part (|101() . requires only the 
classical periods Ti{E) and actions Si{E) of the one-dimensional systems as an input. Nevertheless, 
it contains all information about the periodic orbits of the two-dimensional system - not only the 
degenerate families forming two-dimensional rational tori, but also the existing isolated orbits as 
will be shown explicitly in the following. It also handles all possible bifurcations uniformly. The 
formula (|100|) therefore goes far beyond the standard trace formulae Ej for integrable systems 
which only take the leading rational tori into account and cannot account for bifurcations. 

The integral in (|100)) can, in general, not be done analytically. For a harmonic oscillator 



(102) 



we have Si{E) = 2'KE/uJi and Ti{E) = 2t: /uji {i = u,v). The integral then is elementary and yields 



9ho{E) 



+ 



-iky2iTE /hu!v 



{ky/Uy - ku/<^u) 

For irrational frequency ratios uju ■ ujy, no singularities arise and using the identity |35 



(103) 



sm( z 



E 



(z — kn) 



(z / nvr) 



(104) 



we can sum the first term in (|1()3)) over all ky and the second term over all k^. The result is the 
exact Gutzwiller trace formula for the irrational harmonic oscillator |291 139j which yields its correct 
quantum-mechanical spectrum 
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E 



+ 



1 
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2ttE \ 

huJy J 



(105) 
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The last two terms contain the sums over the (only) isolated periodic orbits along the u and v 
axes. Corresponding trace formulae for rational frequency ratios can be obtained from the above 
by taking suitable limits ■30^. E.g., in the isotropic limit uju = uj^ = uj one obtains 

which is again a quantum-mechanically exact trace formula in terms of the two-fold degenerate 
families of periodic orbits with SU(2) symmetry, having the primitive actions S{E) = 2ttE/lv. Note 
that the standard methods to derive the trace formula for integrable systems HI E] cannot be used 
for harmonic oscillators, since the ("curvature") tensor of second derivatives of the Hamiltonian 
(|1U2|) with respect to the torus actions /j = Ei/coi is identically zero. 

For systems in which the actions Si{E) are no simple functions, the integral in can in 

general only be done numerically. This becomes practically impossible if one wants to generate the 
semiclassical EBK spectrum by summing pOOj) over all ku and k^. In the example treated in the 
next subsection, we show how an asymptotic evaluation of the integral can be used to establish the 
relation to the Berry- Tabor type trace formula for the tori and the Gutzwiller trace formula for 
isolated A orbit, and to derive the same uniform approximation for the bifurcations of the isolated 
A orbit as we have obtained in appendix B using the normal form theory. 



2. Asymptotic evaluation for the separable Henon-Heiles system 
and global uniform approximation 

In the following we specialize to the integrable Henon-Heiles (IHH) system H37() . expressed in the 
scaled coordinates u,v as in (jSHj)- Here = Tb = 27r, Su{E) = Sb{E) = 2TiE, and the period 
Ta{E) and action Sa{E) of the v motion are given in equations - (jSHJ- The TF level density 
is then given by 

9tf{E) = ^Sa[E) (107) 
and the oscillating part can be written as 

5g{E) = ^ y ' / Ta{E^) e^[^-^A(i?.)+2^fc4i?-i?„)]/n ^ ^-^Qg^ 

where the prime indicates that the TF contribution from k^ = k^ = must be left out. Note that 
upon independent summations over k^ and k^ the imaginary parts cancel, consistently with the 
general real expression (|lflO|) . 

In order to evaluate the integral in ()1U8() in the semiclassical limit ?i — > 0, we use the following 
general formula [371 EHl 



rr(x)e^^(-)/^dx ~ ^r(x,) J^^e*^(-»)/^+^Sign[5"(x.)]V4 

Ja . ^ \S {Xi)\ 

+T{h) -Jr- e^5(^)/^— /2 + T{a) e'^ia)/n+in/2 ^ ^^^g-^ 

S'[b) ^ W 

neglecting corrections of higher order in h. Equation (|1()9|) is a generalization of the standard 
stationary-phase integration, taking into account the end-point corrections whose contributions are 
of order h^^'^ relative to those from the stationary points. The sum in the first line above is to be 
taken over all stationary points Xi which lie in the integration interval a < Xi < b. If either of the 
end points a or 6 is a stationary point, its contribution to the sum has to be divided by two and the 
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corresponding term in the second line above has to be omitted. If there is no stationary point at all 
in the interval [a,b], there is no contribution to the first line and the leading terms of the integral 
/ are of order h as given by the end-point contributions in the second line alone. 

The stationarity condition for the phase in the integral (|1()8|) leads to the resonance condition 
for the rational tori 

TAiElkJ = 2vrA:„ = KTb . (110) 

Note that this condition is independent of the energy E. The stationary points = E^^^^ are the 
energies at which the k^-ih. repetition of the A orbit bifurcates, cf. equations (|42|) . (|4.Sj) in section 
4.3. The condition (|ll()j) can only be fulfilled if ku and k^u have the same sign and if j/cul > \k^\- In 
most formulae below, we take ku and k^ to be positive (or zero for one of them) and account for 
the two signs by an extra factor of two in the summations, taking real parts where necessary. For 
ku = kv = k, the stationary point is at £^^^ = 0, for all other tori the stationary points are at finite 
energies. This gives, according to the first line in (|109)) . the following asymptotic contribution to 
lowest order in h: 



oo oo r 1 

1 „ , „s TT 



Sn^k^ (E) + ^ 



e{E-Ei,j, (111) 



which is exactly the Berry- Tabor trace formula Ej with the amplitudes 



_ ku [ 27r ku^ku\ 

-At,„,„ - tKk. ^3/2 y k,T'^{El^,^) ' ^'^-'^ ku = ku]- 

The actions of the tori are 

Sn^.. (E) = ku Sa{EI^O + 2^ku {E - El,J . (E > EIJ (113) 

The diagonal torus T^^, which comes from the lower end point Eu = of the integral in H1U8() . cor- 
responds to the k-th repetition of the classical B orbit. This is somewhat puzzling, since classically 
this orbit appears to be isolated along the u axis, but semiclassically it contributes in the same way 
as the two-dimensional tori with ku ^ ku with an amplitude proportional to fi~^^'^ . The reason for 
this is connected to the fact that the energy E^^f^ = at which it bifurcates is simultaneously the limit 
of the isotropic two-dimensional harmonic oscillator in which all orbits form a two-dimensionally 
degenerate family with SU(2) symmetry. The contributions from fct, = and ku ^ lead to a small 
correction 

*e»'(^)^^£/^s.n(|2.^), (114) 

which is of order y/h with respect to ()112() and found to be negligible in our numerical calculations. 

The upper end point Eu = E of the integral in (|lfl8j) corresponds to motion along the v axis 
which classically gives the isolated A orbit. According to the second line in (|1U9|1 . this yields the 
asymptotic contribution 

( JP) - Ia{EI ^ ^ (-l)fe"+^'^ jlk^SA(E)/h-nm 



Using the identity (|in4|) we can do the summation over ku analytically and find 



s,<.) (E) = ^t '-^'l,, cos sm " 



27r/i ^\u\kuTA{E)l2\ 

Ky i 



h ^ ' 2 



(116) 



which is exactly the Gutzwiller trace formula for the isolated A orbit. Of course, this expression 
cannot be used at the bifurcation energies of the A orbit where, on one hand, H104() is not valid 
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and, on the other hand, the upper end-point correction from the integral (|108)) should be replaced 
by one half of the corresponding torus contribution to The contributions from ku = and 

/c^ 7^ in ()115() lead to a small correction 

*e'(^) = ^E4^s.n(^5.(^,) (117) 

which is included in ()115() . H116|) and will be referred to below. 

We have thus established that the isolated A orbit emerges asymptotically, with its standard 
Gutzwiller amplitude jlUj, from the upper end-point corrections of the EBK trace formula (|1U8() . 
whereas the tori with their standard Berry- Tabor amplitudes come from the stationary points of 
the phase in the integral (|1U8() . 

In order to obtain finite amplitudes at the bifurcations and the symmetry point e = 0, we have 
to develop a uniform approximation. This can be done quite easily by expanding the phase and the 
amplitude of the integrand in ()1U8|) around the stationary points E'^ ^ up to first and second order 
in E'^ — E"^ ^ , respectively. Noting that the torus action St{E) in represents the first two 

terms of the same expansion of Sa{E) around E — E^^^f^^, this leads to the approximate contributions 

5Re|e^^-^.^.(^)-^'^('^"+'=")|^''(6fc„fc„ +Cfe„fc„E,) dE,| (118) 

which are exactly of the same type as those which we have derived from the normal form theory 
in appendix B, and which can be integrated analytically using the formulae H76() and H77() . The 
parameters flfc^fci,, bk^k-u and Cfc„fc„ must be determined by the requirement that the asymptotic 
amplitudes and actions of both the A orbit and the Tkuk^ torus be recovered far away from the 
bifurcation energy E^^^j^^. This procedure is completely analogous to what has been discussed in 
the earlier appendices and need not be repeated here. The final uniform approximation which we 
obtain after summing over all tori is 



oo oo 



kfj 1 ku ky 



h h 



kukv 



cos 



ku kv 



^St, I, {E) 



2 

+ 5g^^'\E) + 5giT\E), (119) 

where we have included the small contributions H114|) and H117() to be consistent up to order in 
the amplitudes, although they are numerically insignificant. Here we have defined 

and 



ik^K = y^^^^ ' A5fc„fc„(i?) = KSa{,E) - St,^,^{E) > 0, (121) 

and the amplitudes and actions of the tori are given in (|112() and (|113() . respectively. 

The uniform trace formula (|119|) is discussed and tested versus quantum-mechanical results in 
section 4.3. Here we compare its results with those of a numerical integration of the EBK trace 
formula H108|) . We choose the value A = 0.04, where the saddle energy is = 104.666, and a 
resolution of the energy spectrum limited by |A;„| < 8. In the figures El -El covering different 
energy regions, the upper panels show the results of (|108() by the solid lines and the sum of the 
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asymptotic Berry- Tabor contributions of the tori (llllj) plus the Gutzwiller contribution (|116|) of 
the isolated A orbit by the dashed lines. In the lower panel, the same results of (|108() (solid lines) 
are compared to those of the uniform approximation H119() (dashed lines). In all cases, the latter 
proves to be an excellent approximation to the exactly integrated EBK trace formula (|l()<Sj) . 
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Figure 16: Shell structure in the level density of the integrable Henon-Heiles system. Solid lines: 
exact EBK integral (|1U8I) . Dashed lines: asymptotic approximations. Upper panel: Berry- Tabor 
contributions of the T tori plus Gutzwiller contribution (|116() of the isolated A orbit. Lower 

panel: uniform approximation (|119|) . In both cases, repetition numbers \ku,k^\ < 8 are included. 
Here the lowest energy region is shown where the isolated A orbit contribution diverges in the limit 
e ^ 0. 



In figure [121 the lowest energy region is shown, where the Gutzwiller contributions of the A orbit 
are seen to diverge in the limit e — > 0. The divergences disappear in the uniform approximation. In 
figure 1171 an intermediate energy region is shown which includes the divergences of the isolated A 
contributions at the bifurcations of the ■ = 9 : 8 and 8:7 resonances. Figure ITHl shows the top 
energy region containing all resonances with ku '■ k^ > 3 : 2. 

These results demonstrate that the uniform approximation (|119|) . which expresses the level 
density in terms of the Berry- Tabor and Gutzwiller amplitudes of the periodic orbits, reproduces 
the numerically integrated EBK trace formula (|lfl8|) to a high degree of accuracy. 

Finally we stress that our above derivation of the uniform approximation H119() is not limited 
to the separable Henon-Heiles system, but can easily be modified to any separable potential by 
starting from the general EBK trace formula H1UU|) rather than from H108|l . 
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Figure 17: Same as figureElin an intermediate energy region. The top panel exhibits the divergence 
of the Gutzwiller contributions of the A orbit (dashed hne) near the : k^j = 9 : 8 and 8:7 
resonances. 
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Figure 18: Same as figure [THl in the top energy range, covering the resonances with k^ : ky > 3 : 2. 
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